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Abstract 

The electromagnetic fields associated with tortuous lightning channels are usually characterised 
by a pronounced fine structure. This work investigates and quantifies the eff'ects of channel tortuos- 
ity on the return stroke electromagnetic field shapes in the close lightning environment (the range 
up to 100m from the lightning striking point). General equations for lightning return stroke elec- 
tromagnetic fields for arbitrarily located observation points are derived from Maxwell's equations. 
In order to include arbitrary channel shapes, the channel is described by a parametric represen- 
tation in Cartesian coordinates with the channel length as the free parameter. The return stroke 
current required for the evaluation of the derived equations is calculated from a current generation 
type model. The field computations show that amplitudes and waveforms of the electromagnetic 
fields in the close lightning environment are considerably infiuenced by the channel shape. In 
particular, the induction component of the electric fields radiated by a tortuous lightning channel 
differs significantly from that associated with a straight channel. 
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I. INTRODUCTION 



In calculating the electromagnetic fields radiated by ground striking lightning discharges, 
the lightning channel is usually idealised by a straight and vertical line between cloud and 
ground. Due to the symmetry in such assumptions, the equations for the electric and 



magnetic fields can be formulated in cylindrical or spherical coordinates [e.g. IsHlOl. Il2l Il3 |. 

The close electromagnetic fields associated with cloud-to-ground lightning, however, is 
characterised by a pronounced fine structure which is due to the tortuosity of the lightning 
channel. Recently, various approaches for the computation of electric and magnetic fields 
radiated by tortuous lightning channels have been published [e.g. 2,0, 16 1. The loss of 



radial symmetry in considering arbitrarily shaped channels necessitates the introduction of 
Cartesian coordinates which implicates an increase of complexity in deriving and solving the 
equations for the electromagnetic fields |16 |. 

Section [TTl of our work presents a detailed derivation of very general expressions for light- 
ning return stroke fields. The equations are formulated in Cartesian coordinates and include 
both arbitrarily shaped lightning channels and arbitrarily located observation points. More- 
over, a height variable return stroke velocity is taken into account. Arbitrary channel shapes 
are considered by the introduction of a parameter representation for the lightning channel 
with the channel length as the free parameter. Such a description allows the direct appli- 
cation of return stroke models which predict the current distribution as a function of the 
length coordinate of the lightning channel. 

In our work, the return stroke current distribution required for the field computations is 
calculated on the basis of the current generation type model proposed by Cooray et al. 
The model includes the attachment process and enables thus to include the upward growing 
connecting leader. Moreover, it allows to take into account a height variable return stroke 
velocity. The model of Cooray et al. ^ and the modifications done in our work are briefly 
discussed in section IIIII 

The equations derived in section Ull are evaluated for both a straight and a tortuous 
lightning channel. The tortuous channel used in the computations is formed from randomly 
generated parameters. The effects of channel tortuosity on the close electromagnetic fields 
are investigated in section UTTI too. 

The presented computations were done within the framework of an investigation on the 
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action of lightning electromagnetic pulses on biological tissue. Biological tissue is mainly 
influenced by the electric flelds which are induced by the transient magnetic pulses associated 
with return strokes. For this reason, a short paragraph of section IIIII deals with the induced 
electric flelds. 



II. DERIVATION OF GENERAL EQUATIONS FOR LIGHTNING RETURN 
STROKE ELECTROMAGNETIC FIELDS 

A. Geometry 

The geometry used in deriving equations for the electric and magnetic flelds radiated 
by a lightning return stroke is shown in flgure [H Cartesian coordinates with basis vectors 
(ej;, e^,, 62) are used. The lightning channel /(s) = lx{s) e^ + ly{s) e^ + /^(s) e^ is parametrised 
by its length s, where s is measured from ground on upward. The base point of the channel is 
assumed to be the point of origin of the coordinate system. R{s) = \J (d/(s) — r) ■ (d/(s) — r) 
denotes the distance between a channel segment d/(s) and the observation point r = xe^; + 
yey + zGz- The velocity distribution of the upward extending return stroke channel is 
assumed to be specifled by a function v{s). Real and retarded channel lengths at time t are 
denoted by s = L{t) and s = Lret(i), respectively. 



B. Solution of Maxwell's equations for the current carrying return stroke channel 

Maxwell's equations in vacuum are given by (since ~ 1, /^r ~ 1 in air, air can be 
idealised as vacuum) 

^.E{f,t) = ^^ (la) 

V-B{f,t) = (lb) 

V X B{f, t) = fioEo + /^o J(r, t) (Id) 

where E{f, t) and B{f, t) denote the electric fleld and the magnetic flux density, respectively. 
For the present case of a current carrying lightning channel parametrised by its length 
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FIG. 1: Geometry used in deriving equations for the electromagnetic fields radiated by a lightning 
return stroke. Real and retarded length of the lightning channel are illustrated by different shades of 
grey. The observation point is marked with f. The grey bar indicates a perfectly conducting ground 
plane. 

coordinate, the current density j{r,t) can be expressed as 

J(r, t) d'r = z{s, t) dtis) = z{s, t) ^ ds (2) 

where i{s,t) denotes the return stroke current distribution. There is no need to replace also 
the charge density p(r, t). This is due to the fact that the scalar potential ${r, t) introduced 
by equation fl3aj) can directly be calculated from the Lorentz gauge (cf. eq. (jlj)). Equations 
(P) are solved by the introduction of the potentials ^(r, t) and A{r,t) defined through 

E(f,t) = -W(f,t)-^^^ (3a) 
5(f, t) = V X A(f, t) (3b) 

Using the Lorentz gauge 
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and combining equations and ([T]), ${r,t) and A{f,t) can be expressed as 

-Lrct(t) t 



<P{r, t) 



1 



ds dr 
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(5b) 



where th = 1^] /c denotes the time taken by a signal to propagate from the hghtning striking 
point to the observation point r, and t^et = t — R{s)/c is the retarded time. Since retardation 
effects have to be taken into account, the s- integration is carried out along the retarded 
channel, that is, the range of integration is given by [s = 0, s = L^ct(t)]. Combining equations 
([3]) and ([5]) yields 
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where the Levi-Civita symbol Skim is given by 



(6b) 



£klr 



1 if {k,l,m) = {x,y,z),{y,z,x),{z,x,y) 
-1 if {k,l,m) = {z,y,x),{y,x,z),{x,z,y) 
otherwise 



(7) 



For the sake of a better oversight the r-, s-, and t-dependences of /(s), R{s), and tret(^, t) 
are omitted in equations ([6]). 



Details to t he p resented solution method of Maxwell's equations can be found in many 



textbooks [e.g. 
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Similar equations for lightning electromagnetic fields are presented in several publications, 
most of them based on the assumptions of a straight and vertical lightning channel and an 
observation point on ground (i.e. z = 0) [e.g. 8- lol. H. 

Zhao and Zhang [16| derived equations for the electromagnetic fields radiated by arbitrar- 
ily shaped lightning channels. Their approach is very similar to that used in our work. In 
contrast to our work, however, their equations do not consider arbitrarily located observation 
points. 

Equations are thus the most general expressions for lightning return stroke fields and 
include all other formulations as special cases. Assuming a straight and vertical lightning 
channel (i.e. l{s) = sBz) and an observation point on ground (i.e. z = 0), for example. 



•g.y, 
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121 1 . On the other 
if an 



equations ([6]) pass into those usually presented in the literature [e.g ^ 

hand, equations ([6]) correspond to the equations presented by Zhao and Zhang [l^ i 
observation point on ground (i.e. 2 = 0) is assumed. 

It should be mentioned that equations ([6]) do not include effects arising from a conducting 
ground plane (cf. subsec. Ill C|) . 



C. The presence of a conducting ground plane 

The presence of earth's surface is not yet included in equations ([6]). Usually, the ground 
plane is taken into account by considering an image channel carrying a current of opposite 
polarity and direction [e.g. llijj. The fields at a given observation point then consist of the 
contributions from the real and the image channel. For an observer on a perfectly conducting 
ground plane the contributions of real and image channel are of equal magnitude, so that the 
vertical electric field and the horizontal magnetic flux density are given by the double of Ez, 
B^, and By from equations ([6]). The horizontal electric field and the vertical magnetic flux 
density vanish on the surface of a perfect conductor. When considering observation points 
above ground, the contributions of real and image channel have to be calculated separately. 
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D. Real and retarded channel lengths for an extending lightning channel 



In order to calculate the real and retarded channel lengths for a return stroke channel 
extending with velocity v{s), the average velocity of the return stroke front between ground 
and length coordinate s is introduced [11 1 



tuis) 



r da 
J v{a) 



(8) 



tu(s) denotes the time taken by the return stroke front to propagate from ground to s. The 
time depending real length L(t) can be calculated by solving the equation 

s=L 

tn{L)= f ^ = t (9) 





Considering retardation effects, that is, taking into account that a given time t comprises 
the time taken by the return stroke front to reach the retarded length Lj-^tit) and the time 
taken by the signal to propagate to the observation point, the retarded channel length can 
be calculated from [cf. lU] 

For observation points above ground the same calculation must be done for the image chan- 
nel, too. 



III. COMPUTATION OF THE CLOSE ELECTROMAGNETIC FIELDS RADI- 
ATED BY STRAIGHT AND TORTUOUS RETURN STROKE CHANNELS 

A. Return stroke model 

Equations ([6]) were evaluated for a current generation type model of the return stroke. The 
model of Cooray et al. [2], which refers to first return strokes and includes the attachment 
process, that is, the preceding upward growing connecting leader, was modified in order to 
describe subsequent return strokes, too. For this purpose, the channel base current used 
by Cooray et al. ^ was replaced by a current waveform typical for subsequent return 
strokes. Moreover, the length of the upward growing connecting leader, which Cooray et 
al. 2! calculated by estimating the potential difference in the gap between connecting and 



dart leader, was replaced by a typical experimental value. In simple terms, the required 
model input parameters are the channel base current, the leader charge distribution, and 
the velocity distribution of both connecting leader and subsequent return stroke. Details 
can be found in Cooray et al. {2]. 

Within the framework of our field computations we investigated also the electromagnetic 
fields radiated by first return strokes. Cooray et al. [2] showed that the close electromag- 
netic fields associated with first return strokes strongly depend on the assumed velocity 
distribution of the return stroke front. The influence of the velocity distribution on the 
electromagnetic fields radiated by subsequent return strokes was not explicitly investigated 
in our work. From a physical point of view, however, the upward growing connecting leader 
can be assumed to have a considerable influence on the close electromagnetic fields, even 
though the connecting leader is comparatively unpronounced for subsequent return strokes. 

The channel base current used in our computations is shown in figure |2l It consists of the 
sum of two Heidler functions [cf. l5| with a parameter set introduced by Diendorfer and Uman 
[sj. The leader charge distribution was assumed to take the constant value of 0.14 mCm~^. 
The length of the connecting leader was set to 10 m. These experimental values were taken 
from Rakov and Uman [6] . The calculated velocity distribution is shown in figure [31 

The return stroke model was implemented in Matlab 7.6 (The MathWorks Inc., Natick, 
US-MA). The analytical formulation and the numerical methods in solving the appearing 



integral equation for the channel current were adapted from Thottappillil and Uman 11| 
and Cooray pL] . The computed return stroke current distribution is shown in figure El 



B. Tortuous lightning channel 

The tortuous lightning channel used in our field computations was formed by a series 
of straight segments. The lengths of the single segments were generated randomly from a 
Gaussian distribution with a mean value of 20 m and a standard deviation of 10 m. The 
inclination angle of each segment, that is, the angle between the segment and the z-axis, 
was taken from a Gaussian distribution with a mean value of 16° and a standard deviation 
of 11.5°. The azimuth was assumed to be uniformly distributed between 0° and 360° . These 
parameters were chosen on the basis of those used by Lupo et al. [4] and Vargas and Torres 



14| . The resulting tortuous channel is shown in figure HI 
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FIG. 2: Return stroke current at different length coordinates in the lightning channel. The waveform 
associated with s = m corresponds to the assumed channel base current. 
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FIG. 3: Return stroke velocity distribution as a function of the length coordinate s. The dashed 
line shows the velocity on a ten times faster length scale where the portion of the upward growing 
connecting leader becomes visible. 

C. Radiated electromagnetic fields 

Return stroke electromagnetic fields were computed for both a straight and vertical, and 
a tortuous channel. In the following discussion the indices "tort" and "str" refer to the fields 
associated with the tortuous and the straight lightning channel, respectively. 

The results for the electric fields and the magnetic flux densities are shown in figures |5] and 
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FIG. 4: 2D and 3D perspectives of the tortuous lightning channel used for the field computations. 
The tortuosity was randomly generated. Note that the lowest channel sections (roughly the first 
50 m) are inclined away from an observer on the positive x-axis and inclined toward an observer 
on the negative x-axis. 

[6l The observation points were assumed to be on ground along the x-axis, that is, |x| gives 
the distance from the hghtning striking point. While for the straight and vertical lightning 
channel the magnetic flux density in x-direction vanishes due to the radial symmetry, the 
tortuous channel radiates small magnetic fields in x-direction, too. Due to the assumption of 
a perfectly conducting ground plane, the horizontal electric fields vanish for both channels. 

According to figures [5] and El the fields associated with the straight and the tortuous 
channel mainly differ in the maximum amplitudes. Moreover, the fields radiated by the 
tortuous channel show an additional fine structure which changes the waveforms of the 
radiated electromagnetic fields. In the following paragraphs these properties are discussed 
in more detail. 

Zhao and Zhang |l6|] presented a similar discussion of the electromagnetic fields radiated 
by tortuous return stroke channels. They investigated the azimuthal dependency of the 
remote fields (100 m — 100 km from the lightning striking point). Our work, on the contrary, 
concentrates on the radial dependency of the very close return stroke fields (20 — 100 m from 
the lightning striking point). Thus, it can be regarded as a supplement to the discussion of 
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Zhao and Zhang [1 



1. Electric fields 

In the hmit t — t- oo the close electric fields shown in figure [5] converge to a constant 
value. This means that they are mainly electrostatic, that is, they are dominated by the 
electrostatic part of equation f pa|) (i.e. the fourth line in eq. f pa|) ). which is the only part 
providing a finite contribution for t — )■ oo. Note that the electrostatic part is proportional to 
R~^, while the other parts of equation f pa|) are proportional to and R~^, respectively. 
Thus, the electrostatic part mainly contributes to the close electric fields. 

The maximum amplitudes of the radiated electric fields considerably depend on the ob- 
servation point (cf. fig. [5]). At a; = 20 m, for example, the maximum amplitude of -Etort 
reaches only 52 % of the value associated with Estr- At x = —20 m, on the other hand, the 
maximum amplitude of -Etort is l.xxx times larger than that of Egti- Both, the increase of 
EtoTt for observation points along the positive x-axis and the decrease of -Etort for observa- 
tion points along the negative x-axis are due to the shape of the tortuous lightning channel: 
The lowest sections (i.e. the first 50 m) of the sample channel shown in figure H] are inclined 
toward an observation point along the positive x-axis and inclined away from an observation 
point on the negative x-axis. After the first 50 m the channel inclination is reversed. This 
is refiected in the electric fields at x = 100 m, where the -Etort is larger than Estr- The close 
electromagnetic fields radiated by tortuous lightning channels thus strongly depend on the 
orientation of those channel segments which primarily contribute to the total fields a given 
observation point. 



2. Magnetic fields 

The orientation of the lowest channel segments is refiected in the magnetic fields, too 
(cf. fig. [6]). At X = 20 m, the maximum value of -Btort reaches roughly 70% of the value 
associated with Bstr- At x = 20 m, on the other hand, the peak value of -Btort is more than 
1.5 times larger than the peak value of -Bgtr. Since the assumed tortuous channel is "quite" 
straight, the main portion (~ 85 %) of the magnetic fields is radiated in ^/-direction (cf. fig. 



6(b), fig. 6(c), and fig. 6(d), fig. 6(e)) 
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FIG. 5: Computed return stroke electric fields for a straight (a) and a tortuous (b) lightning chan- 
nel. The observer was assumed to be located on different positions along the x-axis. Due to the 
assumption of a perfectly conducting ground plane, the horizontal electric fields vanish. 

The magnetic fields associated with the tortuous channel show a pronounced fine struc- 
ture. The relative amplitudes of the additional maxima are observed to increase with in- 



creasing distance from the lightning striking point (most clearly visible in fig. 6(c)). This 
simply refiects the fact that the very close magnetic fields are mainly dominated by the 
lowest channel segments, while with increasing distance also the higher channel sections 
contribute. 

The peak values of -Btort are furthermore observed to be shifted with respect to the peak 



values of -Bstr (cf. fig. 6(a) , fig. 6(b) , fig. 6(d) ). At x = 20 m the peak value of -Btort is reached 
0.3 /is later than that of -Bstr- At x = —20 m it is reached 0.3 /is earlier. The magnitude 
of this shift is observed to decrease with increasing distance. This suggests the shift to be 
due to the orientation of those channel sections which mainly contribute to the fields at a 
given observation point. Generally, a shift of the peak values implicates that -Btort increases 
faster or more slowly than B^tr- If the channel sections contributing to the fields at a given 
observation point are inclined away from the observer, the fields are characterised by a more 
fiat increase. Channel sections inclined towards the observer, on the other hand, result in 
a steeper increase. The channel orientation thus considerably effects the time derivative of 
the magnetic fields. This in turn has a strong infiuence on the induction component of the 
electric fields, that is, the electric fields induced by the transient magnetic pulses. The effect 
of channel tortuosity on the induced electric fields is discussed in subsection IIII C 31 
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FIG. 6: Computed return stroke magnetic flux densities for a straight (a) and a tortuous (h, c, d, e) 
lightning channel. The observer was assumed to be located on different positions along the x-axis. 
The vertical magnetic flux densities of real and image channel cancel on a perfectly conducting 
ground plane. 

3. Induced electric fields 



The results of the field computations presented in the paragraphs IIII CT] and [ITl C 21 were 
used for an investigation on the action of return stroke magnetic pulses on biological tissue. 
Without going into detail it should be mentioned that the action of magnetic fields on bio- 
logical tissue mainly happens by the electric fields which are induced by transient magnetic 
pulses. This paragraph presents a short discussion on the effects of channel tortuosity on 
the induced electric fields. 

The induced electric fields, in the following labelled with E^'^'^, are given by the first 
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Time [us] Time [jis] Time [|j.s] 
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FIG. 7: Computed induced return stroke electric fields for a straight (a) and a tortuous (h, c) 
lightning channel. The observer was assumed to he located on different positions along the x-axis. 
Due to the assumption of a perfectly conducting ground plane, the horizontal induced electric fields 
vanish. 

term of equation fl3al) . Figure [7] shows the results of a numerical evaluation of this term. A 



comparison of figures 7(a) , 7(b) , and 7(c) makes clear that the shape of the lightning channel 



strongly influences the induced electric fields. At x = 20 m the peak value of E^^'^^ reaches 
only 55% of the corresponding value of E^^f. At x = —20 m, on the other hand, the peak 
value of -Etort is more than 1.1% higher than that of E^^f. It is furthermore observed that 
the fields associated with straight and tortuous channels mainly differ in the peak range, 
while their shapes are similar in the initial fast increase and the final slow decrease. 

The induced electric fields reflect the properties of the magnetic fields discussed in sub- 
section nil C 21 The magnetic fields associated with the tortuous channel were found to 
increase faster or more slowly than those associated with the straight channel (remember 
the shifts of the maxima of -Btort)- Since the amplitudes of the induced electric fields depend 
on the time derivative of the magnetic fields (cf. eq. (fTc|) ). this explains the fact that the 
induced electric fields are larger for observation points along the positive x-axis, and smaller 
for observation points along the negative x-axis. The main changes in the induced electric 
fields are due to the fastest changes in the magnetic fields. Therefore, the main changes 
of the induced electric fields concentrate on the first two microseconds. The peak value in 
the induced electric fields corresponds to the point in time with the largest changes in the 
magnetic fields (cf. fig. O fig. [7]). 
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For the sake of completeness it should be mentioned that the induced electric fields 
are of course included in the total electric fields shown in figure [51 In the close lightning 
environment, however, they are superimposed by the electrostatic component and thus not 
explicitly visible. 

D. Computational methods 

Equations (E]) were implemented in Matlab 7.6 (The MathWorks Inc., Natick, US-MA). 
In carrying out the numerical integrations the trapezoidal rule was used. For the field 
integrations a total channel length 5* (cf. fig. [1]) of 1000 m was assumed, which is enough 
when the fields on ground within 100 m from the lightning striking point are considered. 
Time and length steps were set to 5 ■ 10"^ s and 1 m, respectively. For the fields at distances 
of 20 m and 40 m partly also higher resolutions were needed. 

IV. SUMMARY AND CONCLUSIONS 

Very general equations for lightning electric and magnetic fields including both arbitrar- 
ily shaped lightning channels and arbitrarily located observation points were derived from 
Maxwell's equations. The lightning channel was described by a parametric representation in 
Cartesian coordinates. Equations for the determination of real and retarded channel lengths 
for the case of a height variable return stroke speed were presented. In order to discuss 
the effect of channel tortuosity on the close return stroke electromagnetic fields, the derived 
equations were solved for a current generation type model of the lightning return stroke. 
Return stroke electromagnetic fields were computed for both a straight and vertical, and a 
randomly generated tortuous lightning channel. Moreover, the electric fields induced by the 
transient return stroke magnetic pulses were investigated. 

The results show that the global amplitudes of the close electromagnetic fields radiated by 
a tortuous lightning channel mainly depend on the orientation of the lowest channel sections. 
The relative amplitudes of the additional fine structure due to the channel tortuosity are 
observed to increase with increasing distance from the striking point, which was reduced to 
the fact that with increasing distance more different orientated channel sections contribute 
to the fields at a given observation point. Furthermore, the maximum of the magnetic 
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fields radiated by tortuous lightning channels was found to be shifted with respect to the 
magnetic fields associated with a straight channel. The shift implicates a change in the time 
derivative of the magnetic fields radiated by tortuous channels. This in turn was found to 

have a remarkable influence on amplitudes and waveforms of the induced electric fields. 
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